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On the basis recent seven-loop perturbation expansion for u^ 1 — 3/(2 — a) we perform a careful 
reinvestigation of the critical exponent a governing the power behavior \T C — T\~ a of the specific 
heat of superfluid helium near the phase transition. With the help of variational strong-coupling 
theory, we find a = —0.01126 ±0.0010, in very good agreement with the space shuttle experimental 
value a = -0.01056 ± 0.00038. 
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1. The critical exponent a characterizing the power behavior \T C — T\~ a of the specific heat of superfluid helium near 
the transition temperature T c is presently the best-measured critical exponent of all. A microgravity experiment in 
the Space Shuttle in October 1992 rendered a value with amazing precision p) 
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a ss = -0.01056 ±0.00038. (1) 



This represents a considerable change and improvement of the experimental number found a long time ago on earth 
byG. Ahlers §: 

m ! 

>. a = -0.026 ± 0.004, (2) 

r->; 

in which the sharp peak of the specific heat was broadened to 10 6 K by the tiny pressure difference between top and 
bottom of the sample. In space, the temperature could be brought to within 10 _8 K close to T c without seeing this 
broadening. 

The exponent a is extremely sensitive to the precise value of the critical exponent v which determines the growth 
of the coherence length when approaching the critical temperature, £ ex \T — T c \~ v . Since v lies very close to 2/3, and 
a is related to v by the scaling relation a = 2 — 3v, a tiny change of v produces a large relative change of a. Ahlers' 
C^ ' value was for many years an embarrassment to quantum field theorists who never could find a quite as negative — 
the field theoretic z^-value came usually out smaller than y^hl = 0.6753 ±0.0013. The space shuttle measurement was 
1 ■ therefore extremely welcome since it comes much closer to previous theoretical values. In fact, it turned out to agree 
extremely well with the most recent theoretical determination of a by strong-coupling perturbation theory [Bj based 
on the recent seven- loop power series expansions of v (|,|7],|[ , which gave Q 



a sc = -0.0129 ±0.0006. (3) 



The purpose of this Letter is to present yet another resummation of the perturbation expansion for v~ l and for 
a = 2 — iv by variational perturbation theory, applied in a different way than in (q| . Since it is a priori unclear which 
of the two results should be more accurate, we combine them to the slightly less negative average value with a larger 
error 

a sc = -0.01126 ±0.0010. (4) 

Before entering the more technical part of the paper, a few comments are necessary on the reliability of error 
estimates for any theoretical result of this kind. They can certainly be trusted no more than the experimental 
numbers. Great care went into the analysis of Ahlers' data 0. Still, his final result (||) does not accommodate the 
space shuttle value (Fy). The same surprise may happen to theoretical results and their error limits in papers on 
resummation of divergent perturbation expansions, since there exists so far no safe way of determining the errors. 
The expansions in powers of the coupling constant g are strongly divergent, and one knows accurately only the first 
seven coefficients, plus the leading growth behavior for large orders k like r y(—a) k k\kT(k + b). The parameter b is 
determined by the number of zero modes in a solution to a classical field equation, a is the inverse energy of this 
solution, and 7 the entropy of its small oscillations. 
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The shortness of the available expansions and their divergence make estimates of the error range of the result a 
rather subjective procedure. All publications resumming critical exponents such as a calculate some sequences of 
./Vth-order resummed approximations cxn, and estimate an error range from the way these tend to their limiting 
value. While these estimates may be statistically significant, there are unknown systematic errors. Otherwise one 
should be able to take the expansion for any function /(<?) = f(a(g)) and find a limiting number f(a) which lies in 
the corresponding range of values. This is unfortunately not true in general. Such reexpansions can approach their 
limiting values in many different ways, and it is not clear which yields the most reliable result. One must therefore 
seek as much additional information on the series as possible. 

One such additional information becomes available by resumming the expansions in powers of the bare coupling 
constant go rather than the renormalized one g. The reason is that any function of the bare coupling constant /(go) 
which has a finite critical limit approaches this limit with a nonlcading inverse power of g^ , where uj is called the critical 
exponent of approach to scaling, whose size is known to be about 0.8 for supcrfluid helium. Any resummation method 
which naturally incorporates his power behavior should converge faster than those which ignore it. This incorporation 
is precisely the virtue of variational perturbation theory, which we have therefore chosen for the resummation of a. 

For a second additional information we take advantage of our theoretical knowledge on the general form of the 
large-order behavior of the expansion coefficients 



7 (-a) fc fc!fcr(fc + 6) M + 
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In the previous paper pi we have done so by choosing the nonleading parameters d to reproduce exactly the first 
seven known expansion coefficients of a. The resulting expression (0) determines all expansion coefficients. The so- 
determined expression (pj) predicts approximately all expansion coefficients, with increasing precision for increasing 
orders. The extended power series has then been resummed for increasing orders N, and from the ./V-behavior we 
have found the a- value ^) with quite a small error range. 

As a third additional information we use the fact that we know from theory |3j in which way the infinite-order 
result is approached. Thus we may fit the approximate values ajv by an appropriate expansion in 1/iV and achieve 
in this way a more accurate estimate of the limiting value than without such an extrapolation. The error can thus 
be made much smaller than the distance between the last two approximations, as has been verified in many model 
studies of divergent series [|| . 

The strategy of this paper goes as follows: We want to use all the additional informations on the expansion of the 
critical exponent a as above, but apply the variational resummation method in two more alternative ways. First, 
we reexpand the series a(go) in powers of a variable h whose critical limit is no longer infinity but h = 1. The 
closer distance to the expansion point h = leads us to expect a faster convergence. Second we resum two different 
expansions, one for a, and one for f(a) = u^ 1 = 3/(2 — a). From the difference in the resulting a-values and a 
comparison with the earlier result (ra) we obtain an estimate of the systematic errors which is specified in Eq. (0) . 

2. The seven-loop power series expansion for v in powers of the unrenormalizcd coupling constant of 0(2)-invariant 
<^ 4 -theory which lies in the universality class of supcrfluid helium reads [|4| 

v- 1 = 2 - 0.4.go + 0.4681481481482289 g 2 - 0.66739 g 3 + 1.079261838589703 g Q 4 - 1.91274 g 5 

+ 3.644347291527398 g 6 - 7.37808 g a 7 + . . . . (6) 

By fitting the expansion coefficients with the theoretical large-order behavior (0), this series has been extended to 
higher orders as follows g 

Az^ 1 = 15.75313406543747 go 8 - 35.2944 ff0 9 + 82.6900901520064 g a 10 - 202.094 5o n + 514.3394395526179 g 12 

- 1361.42 5o 13 + 3744.242656157152 g 14 - 10691.7 5o 15 + • • • • (J) 

The renormalized coupling constant is related to the unrenormalized one by an expansion g = X)fc=i a feS'o- Its 
power behavior for large go is determined by a series 

s = °, g0 ^ = 1 - .go + 947g " - 2.322324349407407 ff0 3 + 4.276203609026057 g 4 
a log go 675 

-8.51611440473227 ff o 5 + 18.05897631325589 ffo 6 + --- ■ (8) 

A similar best fit of these by the theoretical large-order behavior extends this series by 
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As = 40.38657228730114 5o 7 + 94.6453399123477 g s - 231.3922442162566 g 9 + 588.3206172579102 g 
- 1552.116358404217 g n +4242.372685080157 g 12 - 12001.18866491822 g 13 + 35115.23006646194 
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106234.4643086436 g 15 + 332239.2175082959 g 16 



(9) 



Scaling implies that 3(170) becomes a constant for go — > 00, implying that the power s goes to zero in this limit. By 
inverting the expansion for s, we obtain an expansion for v~ x in powers of h = 1 — s as follows: 

z/ _1 (/i) = 2-0.4 fr-0.093037 fr 2 +0.000485012 /i 3 -0.0139286 /i 4 +0.007349 /i 5 -0.0140478 fr 6 +0.0159545 fr 7 -0.029175 h s 

+ 0.0521537 /i 9 -0.102226/i 10 + 0.224026/i n -0.491045/i 12 + l. 22506 h 13 -3. 00608 /i 14 +8.29528/i 15 -22. 5967 h 16 . (10) 

This series has to be evaluated at h = 1. For estimating the systematic errors of our resummation, we also calculate 
from (III]) a series for a — 2 — 2>v 

a(h) = 0.5-0.3 /i-0.129778/i 2 -0.0395474/i 3 -0.0243203/i 4 -0.0032498/i 5 -0.0121091/i 6 

+ 0.00749308 h 7 -0.0194876 /i 8 +0.0320172 /i 9 -0.0651726 h 10 +0.14422 /i n -0.315055 h 12 

+ 0.802395/j 13 -1.95455/i 14 +5.49143/i 15 -14.8771/i 16 + ... . (11) 

3. In order to get a rough idea about the behavior of the reexpansions in powers of h, we plot their partial sums at 
h = 1 in the upper row of Fig. [j] 
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FIG. 1. Upper plots: Results of partial sums of series (FEfl) for v^ 1 up to order N, once plotted as un = 1/v^ , and once as 
cw = 2 — 3vn- The third plot shows the corresponding partial sums of the series for a. The dotted line is the experimental 
space shuttle value a BB of Eq. (hi). Lower plots: The corresponding resummed values and a fit of them by cq + ci/N 2 + C2/N 4 . 
The constant Co is written on top, together with the seventh-order approximation (in parentheses). The square brackets on top 
of the left-hand plot for v shows the corresponding a- values. 

After an initial apparent convergence, these show the typical divergence of perturbation expansions. 
A rough resummation is possible using Pade approximants. The results are shown in Table |. The highest Pade 
approximants yield 
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a raa = -0.0123 ±0.0050. 
The error is estimated by the distance to the next lower approximation. 



(12) 



4. We now resum the expansions v (h) and a(h) by variational perturbation theory. This is applicable to divergent 
perturbation expansions 



/(*) = X! a » a 



(13) 



which behave for large x like 



f(x) = a*/' ^ 6r, 



-2m/ q 



(14) 



It is easy to adapt our function to this general behavior. Plotting the successive truncated power series for v 1 (h) 
against h in Fig. |2j, we see that this function will have a zero somewhere above h = ho = 3. 
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FIG. 2. Successive truncated expansions of f -1 (/i) of orders N — 2, . . . , 12. 

We therefo re g o over to the variable x defined by h = h(x) = hox/(ho — 1 + x), in terms of which f(x) = v~ 1 {h{x)) 
behaves like (0) with p = and q = 2, and has to be evaluated at aj = 1. The large- a; behavior is imposed upon 
the function with the expansion (O) as follows. We insert an auxiliary scale parameter k and define the truncated 
functions 



N 
f N (x)=K,Pj2 a n(^) ■ 



(15) 



n=0 



The parameter k will be set equal to 1 at the end. Then we introduce a variational parameter K by the replacement 

K -> VK 2 + k 2 - K 2 . (16) 

The functions Jn(x) are so far independent of K . This is changed by expanding the square root in (|16|) in powers of 



K , thereby treating this difference as a quantity of order x. This transforms the terms n p x n / ' n qn in (|15[) into 



polynomials of r = (k — K 2 )/K 



K qn 



k p 



Km 



(p - ?n)/2 



m)/2\ f{p-qn)/2\ f{p-qn)/2\ N „ r 

1 ) r+ { 2 J r +---+I, AT-n ) r 



(17) 



Setting now k = 1, and replacing the variational parameter K by v defined by K 2 = x/w, we obtain from (15) at 
x = 1 the variational expansions 
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f; a„^-^ 2 [i + (« - i)]£:r )/2 , 

n=0 



(18) 



where the symbol [1 + A\^_ q ™ is a short notation for the binomial expansion of (1 + A)^ p «™)/ 2 in powers of A up 
to the order A N ~ n . 

The variational expansions are optimized in v by minima for odd, and by turning points for even N, as shown in 
Fig. ||. The extrema are plotted as a function of the order N in the lower row of Fig. [I[ The left-hand plot shows 
directly the extremal values of v~^ l {v), the middle plot shows the a- values ajv = 2 — 3z/jv corresponding to these. The 
right-hand plot, finally, shows the extremal values of oin(v). All three sequences of approximations are fitted very well 
by a large N expansion c + ci/N 2 + c 2 /A f4 , if we omit the lowest five data points which are not yet very regular. The 
inverse powers 2 and 4 of N in this fit are determined by starting from a more general ansatz cq + c\/N Pl + C2/N P2 
and varying p\,p2 until the sum of the square deviations of the fit from the points is minimal. 
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FIG. 3. Successive variational functions v^Qi) and ajvCi) with TV = 3, . . . , 12 of Table O plotted for /i = x = 1 against 
the variational parameter 7^ = Wx/v, together with their minima for odd N, or turning points for even N. These points are 
plotted against TV" in the lower row of Fig. hi, where they are extrapolated to TV — > oo, yielding the critical exponents. 

The highest-order data point is taken to be the one with N = 12 since, up to this order, the successive asymptotic 
values Co change monotonously by decreasing amounts. Starting with N — 13, the changes increase and reverse 
direction. In addition, the mean square deviations of the fits increasing drastically, indicating a decreasing usefulness 
of the extrapolated expansion coefficients in (R) and (Bh for the extrapolation N — > oo. From the parameter Co of 
the best fit for a which is indicated on top of the lower right-hand plot in Fig. |l|, we find the critical exponent 
a = —0.01126 stated in Eq. (H), where the error estimate takes into account the basic systematic errors indicated by 
the difference between the resummation of a — 2 — 3i^, and of v~ l , which by the lower middle plot in Fig. |l] yields 
a = —0.01226. It also accommodates our earlier seven-loop strong-coupling result (H) of Ref . H. The dependence on 
the choice of h is negligible as long as the resummed series v~ l (x) and a(x) do not change their Borel character. 
Thus ho = 2.2 leads to results well within the error limits in (Eh. 

Our number as well as many earlier results are displayed in Fig. 0. The entire subject is discussed in detail in the 
textbook H. Kleinert and V. Schulte-Frohlinde, Critical Exponents from, Five-Loop Strong- Coupling (j) 4 -Theory in J^-e 
Dimensions , World Scientific, Singapore, 2000 ( http://www.physik.fu-bcriin.de/ ~klcincrt/rc.html=#=bg| ) 
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Note added in proof: 

A recent calculation of a by an improved high-temperature expansion yields the exponent a = —0.0150(17) [M. 

Campostrini, A. Pelissetto, P. Rossi, and E. Vicari, Phys. Rev. B 61, 5905 (2000)]. 
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FIG. 4. Survey of experimental and theoretical values for a. The latter come from resummed perturbation expansions of 
(//-theory in 4 — e dimensions, in three dimensions, and from high-temperature expansions of XY-models on a lattice. The 
sources are indicated below. 
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TABLE I. Results of the Pade approximations Pmat(/i) at h — 1 to the power series v 1 (h) and a(h). The parentheses show 
the associated values of a and v. 
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v (a) 


(u) a 


44 
5 4 

4 5 

5 5 


0.678793 (-0.0363802) 
0.671104 (-0.0133107) 
0.670965 (-0.0128940) 
0.670756 (-0.0122678) 


(0.678793) -0.0363802 
(0.670965) -0.0128940 
(0.670901) -0.0127031 
(0.670756) -0.0122678 



TABLE II. Variational reexpansions of vjf (h) and ajv(ft) for N = 2, . . . ,9 at h — x — 1 which are plotted in Fig. tfl and 
whose minima and turning points are extrapolated to N == oo in the lower left- and right-hand plots of Fig. hi. The lists are 
carried only to N — 9, to save space, whereas the plots are for N — 3, ... ,12. 
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